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ABSTRACT 

We numerically investigate the dynamical and chemical processes of the formation of elliptical galaxies 
in a cold dark matter (CDM) universe, in order to understand the origin of the mass-dependence of 
the photometric properties of elliptical galaxies. Our three-dimensional TREE N-body/SPH numerical 
simulations of elliptical galaxy formation take into account both Type II (SNe II) and Type la (SNe la) 
supernovae (SNe) and follow the time evolution of the abundances of several chemical elements (C, O, Ne, 
Mg, Si, and Fe). Moreover we compare different strengths of SNe feedback. In combination with stellar 
population synthesis, we derive the photometric properties of simulation end-products, including the 
magnitude, color, half-light radius, and abundance ratios, and compare them with the observed scaling 
relations directly and quantitatively. We find that the extremely strong influence of SNe is required 
to reproduce the observed color-magnitude relation (CMR), where we assume each SN yields energy of 
4 x 10 51 ergs and that 90% of this energy is ejected as kinetic feedback. The feedback affects the evolution 
of lower mass systems more strongly and induces the galactic wind by which a larger fraction of gas is 
blown out in a lower mass system. Finally higher mass systems become more metal rich and have redder 
colors than lower mass systems. We emphasize based on our simulation results that the galactic wind is 
triggered mainly by SNe la rather than SNe II. In addition we examined the Kormendy relation, which 
prescribes the size of elliptical galaxies, and the [Mg/Fe]-magnitude relation, which provides a strong 
constraint on the star formation history. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: formation — galaxies: evolution 
— galaxies: stellar content 



1. INTRODUCTION 

Despite the large variety in the internal dynamics and 
structures, elliptical galaxies follow the various scaling re- 
lations. These relationships are expected to contain valu- 
able information about physical processes of their forma- 
tion. For example, colors of elliptical galaxies strongly 
correlate with their luminosities (e.g., Sandage 1972; 
Aaronson, Persson, & Frogel 1981; Bower, Lucey, &; El- 
lis 1992a,b, hereafter BLE92a and BLE92b). This rela- 
tion is usually called the color-magnitude relation (CMR). 
The CMR is conventionally interpreted as an effect of the 
galactic metallicity increasing with the mass (e.g., Faber 
1973; Dressier 1984), although the age difference also con- 
tributes on the CMR possibly (Worthey, Trager, & Faber 
1996). This interpretation is supported by line index anal- 
ysis of galaxies (Kuntschner & Davies 1998; Kuntschncr 
1998; Terlevich et al. 1999; Kuntschner 2000) and slopes 
of CMRs observed in the intermediate redshift (z ~ 0.8) 
clusters (Kodama & Arimoto 1997). The mass-metallicity 
sequence can be naturally explained by the galactic wind 
scenario (Larson 1974) and reproduces the CMR very 
well, according to the previous studies employing the evo- 
lutionary method of population synthesis (e.g., Arimoto & 
Yoshii 1987; Matteucci & Tornambe 1987; Gibson 1997; 
Tantalo et al. 1998). In this scenario, galactic wind is 
induced progressively later in more massive galaxies ow- 
ing to deeper potential and stellar populations in brighter 
galaxies are much more enhanced in heavy elements, lead- 



ing to redder colors. However, these studies assumed a 
simple galaxy model which ignored the internal structure 
and complex star formation history in a forming galaxy. 

Recently developed semi-analytic model (White & Rees 
1987; White & Frenk 1991) makes it possible to describe 
some complex star formation histories, such as a star-burst 
induced by a coalescence of galaxies. This model is based 
on the hierarchical clustering scenario expected in the cold 
dark matter (CDM) cosmology. The formation and merg- 
ing histories of dark matter halos are deduced by an ex- 
tended Press-Schechter theory (Bower 1991; Bond et al. 
1991) or a cosmological N-body simulation. The model 
presumes that gas initially cools into a disk and forms 
stars inside the halo. A spheroidal stellar system is as- 
sumed to be formed by a major merger of two disk galaxies 
(Kauffmann, White, & Guiderdoni 1993; Baugh, Cole, 
& Frenk 1996). Kauffmann & Chariot (1998) showed 
by taking account of chemical enrichment that the CMR 
of elliptical galaxies is reproduced in their semi-analytic 
model with the assumption that the supernova feedback is 
much efficient. Although the recent semi- analytic models 
seem to succeed in explaining the observational properties 
of galaxies at various redshifts, they inevitably assumed 
a phcnomenological model involving a number of parame- 
ters to describe the processes of radiative cooling, star for- 
mation and feedback within the galactic halo. Moreover, 
structures of the end-products could not be discussed, be- 
cause there is no information about the dynamical history 
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of gas and star components. 

Numerical simulations are powerful tools to treat the 
complex physical processes in galaxy formation self- 
consistently. In recent works using three dimensional nu- 
merical simulations which calculate the dynamical and 
chemical evolution self-consistcntly, photometric proper- 
ties of end-products were investigated in combination with 
the stellar population synthesis and compared to the ob- 
servations directly (e.g., Bekki & Shioya 1998, 1999; Con- 
tardo, Steinmetz, & Fritze-von Alvensleben 1998; Mori, 
Yoshii, & Nomoto 1999; Steinmetz & Navarro 1999; 
Koda, Sofue, & Wada 2000a, b; Navarro & Steinmetz 
2000; Kawata 1999, 2001, hereafter K99, K01). In order 
to understand the connection between the observed scal- 
ing relations and the physical processes of the formation of 
elliptical galaxies, we investigate the dynamical and chem- 
ical evolution of elliptical galaxies in a cold dark matter 
(CDM) universe using the numerical simulation. In this 
paper we focus on understanding the origin of the CMR, 
i.e., how and what physical processes control the mass de- 
pendence of the photometric properties. To this end, we 
adopt the elliptical galaxy formation model of K99. This 
model is based on a semi-cosmological model proposed by 
Katz & Gunn (1991) and treats a collapse of a top-hat 
over-dense sphere which initially follows a Hubble flow ex- 
pansion and has a solid-body rotation as an effect of the 
external tidal field. In addition this sphere includes small- 
scale density fluctuations expected in a CDM universe. 
The galaxy is built by merger of small clumps formed due 
to the small-scale fluctuations, i.e., a hierarchical cluster- 
ing. K99 studied the evolution of a seed galaxy which 
has a slow rotation corresponding to a spin-parameter of 
A = 0.02 and showed that the end-product reproduces the 
observed properties of elliptical galaxies very well. Cos- 
mologically, a small spin-parameter is preferred in a galac- 
tic halo collapsing at a high redshift (Heavens & Peacock 
1988). The reason for this is that higher peaks have a 
shorter collapse time; thus, they have less time to get 
spun up by the environmental tidal force. Moreover, a 
significant fraction of elliptical galaxies are considered to 
be the systems which had collapsed at high rcdshifts, ac- 
cording to the zero-point and tightness of the CMR at 
the intermediate redshift (e.g., Bower, Kodama, & Ter- 
levich 1998; van Dokkum et al. 2000). It is thus a nat- 
ural supposition that some elliptical galaxies were formed 
in a halo collapsing at a high redshift and rotating slowly 
(e.g., Blumcnthal et al. 1984). Although K99 studied the 
dynamical and chemical evolution of a galaxy with a fixed 
mass, we here examine the photometric properties of the 
end-products with different masses. We pay a particu- 
lar attention to the supernovae feedback, which is one of 
the most difficult processes to model in the present-day 
numerical simulations, but considered to be a crucial de- 
terminant of the nature of the stellar population of the 
end-products. Therefore, we examine how the difference 
in the strength of the feedback affects the final CMR. In 
addition our numerical code takes account of the chem- 
ical and dynamical feedback of both Type II (SNe II) 
and Type la (SNe la) supernovae (SNe) and follows the 
evolution of the abundances of several chemical elements 
(C, O, Ne, Mg, Si, and Fe). We evaluate relative impor- 
tance of two types of SNe, SNe II and SNe la, on the 



dynamical and chemical evolution. Since our numerical 
simulations naturally provide the information about the 
structure of end-products with enough spatial resolution, 
we compare simulation results not only with the observed 
CMR but also with the Kormendy relation, which pre- 
scribes the size of elliptical galaxies as a function of their 
luminosity. Due to self-consistent calculation of evolution 
of chemical elements abundances, we can also study the 
mass dependence of abundance ratios. Kuntschner (1998, 
2000) studied early-type galaxies in the Fornax cluster, us- 
ing the line indices analysis for high-quality data. He found 
that Fornax ellipticals are mainly old and also discovered 
a strong relation between [Mg/Fe] and the velocity dis- 
persion (see also Trager et al. 2000b). J0rgcnscn (1999) 
studied Coma elliptical and obtained a strong correlation 
between [Mg/Fc] and not only the velocity dispersion but 
also the luminosity. Hereafter we call the correlation be- 
tween [Mg/Fe] and the luminosity "[Mg/Fc]-magnitudc 
relation". Because Mg and Fe are mostly produced by 
SNe II and SNe la respectively and because SNe la have 
a longer delay than SNe II after formation of stars, these 
correlations give strong constraints on the star formation 
history of elliptical galaxies (Matteucci, Ponzone, & Gib- 
son 1998; Thomas, Greggio, & Bender 1999). We ex- 
amine the [Mg/Fe]~magnitude relation for the simulation 
end-products and compare it with the observation. 

The plan for the remainder of this paper is as follows. 
The numerical method and the model used in this paper 
are described in Sections 2 and 3. We show the procedure 
to derive the photometric properties of the simulation end- 
products in Section 4. In Section 5, results of numerical 
simulations are presented and are compared with the ob- 
served global scaling relations. Our discussion and conclu- 
sions are given in Section 6. 

2. THE CODE 

Numerical simulations are performed using an update 
version of the code in K99 so that SNe la feedback is 
taken into account. Because we have already described 
the details of the code in K99, here we describe mainly 
the modelings of star formation and their feedback which 
are revised. Our code is essentially based on the TreeSPH 
(Hernquist & Katz 1989; Katz, Weinberg, & Hernquist 
1996), which combines the tree algorithm (Barnes & Hutt 
1986) for the computation of the gravitational forces with 
the smoothed particle hydrodynamics (SPH: Lucy 1977; 
Gingold & Mnoraghan 1977) approach to numerical hy- 
drodynamics. The dynamics of the dark matter and stars 
is calculated by the N-body scheme, and the gas compo- 
nent is modeled using the SPH. It is fully Lagrangian, 
three-dimensional, and highly adaptive in space and time 
owing to individual smoothing lengths and individual time 
steps. Moreover, it self-consistently includes almost all the 
important physical processes in galaxy formation, such as 
self-gravity, hydrodynamics, radiative cooling, star forma- 
tion, supernova feedback, and metal enrichment. 

The radiative cooling which depends on the metallicity 
(Theis, Burkert, & Hensler 1992) is taken into account. 
The cooling rate of a gas with the solar metallicity is larger 
than that for a gas of the primordial composition by more 
than an order of magnitude. Thus, the cooling by metals 
should not be ignored in numerical simulations of galaxy 
formation (Kaellander & Hultman 1998; Kay et al. 2000). 
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2.1. Star Formation 

We model star formation using a method similar to 
that of Katz (1992) and Katz, Weinberg, & Hcrnquist 
(1996). We use the following three criteria for star forma- 
tion: 1) the gas density is greater than a critical density, 
Pcrit = 2 x 10~ 25 g cm~ 3 , i.e., ~ 0.1cm -3 , following 
Katz, Weinberg, & Hcrnquist (1996); 2) the gas veloc- 
ity field is convergent, V • Vi < 0; and 3) the Jeans un- 
stable condition, h/c s > t g , is satisfied, here h, c s , and 
t g = ^/3-7r/16Gp g are the SPH smoothing length, the 
sound speed, and the dynamical time of the gas respec- 
tively (see K99 for details). The Jeans condition was ig- 
nored in K99 and KOI, because the gas whose density is 
higher than /9 cr ; t usually satisfies this condition. Since the 
Jeans unstable condition might be important in the low 
mass systems, we take it into account in this paper. 

When a gas particle is eligible to form stars, its star- 
formation rate (SFR) is 

dp* = _df>s = c *Pg 
dt dt t g ' 

where c* is a dimensionless SFR parameter and t s is the 
dynamical time, which is longer than the cooling timescale 
in the region eligible to form stars. This formula corre- 
sponds to the Schmidt law that SFR is proportional to 
Pg 5 . The value of c* controls the SFR with respect to the 
local gas density, namely star formation efficiency. When a 
large c* is assumed, many stars can be formed well before 
the system collapses completely, so that the effective ra- 
dius of the final stellar system becomes large. In this paper 
we assume for simplicity that is constant irrespective of 
the mass of the system and calibrated c* so that model 
A2 defined in Section 3 could reproduce the observed Ko- 
rmendy relation. Figure 11 shows the results of model A2 
in the case of c* = 1 (cross) and c* = 0.5 (triangle just un- 
der the cross) respectively. We can see that c* = 1 leads to 
a little larger effective radius than the observed elliptical 
galaxies under the assumed distance modulus (see Section 
5 for details) . Although the difference in c* a little affects 
the other photometric properties which we examine in this 
paper (Fig. 5 and 12), finally we found that this effect is 
quite small, compared to difference in the strength of SNc 
feedback which we focus on in this paper. Consequently 
we set c* = 0.5. In addition the above critical gas den- 
sity, pcrit, is higher than that in K99 and KOI. We also 
carried out a simulation with the same p cr i t as K99 and 
found that this difference does not change the effective ra- 
dius in model A2, because the gas density becomes higher 
than these densities very quickly once the system collapses. 
The difference in p cr it within this range does not affect our 
results. 

We assume that the stars which are represented by a star 
particle are distributed according to the Salpeter (1955) 
initial mass function (IMF). The IMF by number, $(m), 
in each mass interval dm is defined as 

$(m)rfm = Am- (1+x) dm, (2) 

where x = 1.35 is the Salpeter index (this corresponds to 
an IMF by mass cx m~ x ) and the coefficient A is deter- 
mined by the normalization in the mass range M\ < m < 
M u . We set M u = 60 M Q and M\ = 0.4 M Q respectively. 

2.2. The Feedback 



We take account of the energy feedback and metal en- 
richment to the surrounding gas by SNe. We consider here 
both SNc II and SNe la. The simulation follows the evo- 
lution of the abundances of several chemical elements (C, 
O, Ne, Mg, Si, and Fe). 

2.2.1. Type II Supernovae 

For simplification, we assume that each massive star 
(> 10 M ) explodes as a type II supernova within the 
simulation time step in which it was born (instantaneous 
recycling). From equation (2) the number of SNe II in 
a new-born star particle of m s Mq becomes iVsNcii = 
9.26 x 10 _3 m s . The total stellar mass contributing to 
SNe II is M SNo ii = 1.83 x 10~ 1 m s M and the mass of 
(MsncII — ITAsNoii) M Q is returned to the the surround- 
ing gas, where the remnant mass is assumed to be 1.4 M . 
In this process, a portion of the metals produced in the 
stars is also returned to the gas, leading to chemical en- 
richment. We use the stellar yields derived by Notomo et 
al. (1997a) who provide for SNe II the z-th heavy element 
mass, M i; sNeii("i), produced in a star of main-sequence 
mass of m. Then the synthesized mass of SNe II in a new- 
born star particle of m s M for each element is calculated 

by 

M^sNdi = m s Mi^uim^mr^'^dm / / m' x dm M 

J 10 M / J Mi 

(3) 

Table 1 shows M^sNeii for each element. Here we assume 
that the heavy-element production from a star of 10 M 
is negligible, following Tsujimoto et al. (1995). 

These mass, energy, and heavy elements are smoothed 
over the neighboring gas particles using the SPH smooth- 
ing algorithm. For example, when the i-th particle changes 
from gas to a star, the increment of the mass of the j'-th 
neighbor particle due to explosion of the new-born star is 
given by 

77% ■ 

AMsnj = — M SNti W(rij, hij), (4) 

Pg,i 

where 

P e ,i = (p K (xi)) = m j W(r ij J> ij ) (5) 

and W(rij, hij) is an SPH kernel (see K99). 

2.2.2. Type la Supernovae 

Comparing the predicted nucleosynthesis products of 
SNc la and SNe II with the solar abundances of heavy 
elements and their isotopes, Tsujimoto et al. (1995) de- 
rived that the ratio of the total numbers of SNe la to SNe 
II is NsNeia/NsNeii = 0.15. If the progenitor of SNe 
la is a binary system one of which has a mass range of 
m = 3-8 M , AfsNoia/^VsNoii is calculated by 

Nsmu Jwu° m - {1+X)dm 

Thus, we adopt A\ a ~ 0.04. Then the number of SNe la 
in a new-born star particle of m s M becomes AgNoia = 
1.52 x 10~ 3 m s and the ejected mass is assumed to be A\ a 
times the total mass of stars whose mass is between 3 M 
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to 8 M©, i.e., MsNda = 6.94 x 1CT 3 to s M . The nucle- 
osynthesis prescriptions for SNe la are taken from Notomo 
et al. (1997b) as shown in Table 1. If SNe la occur in a 
binary systems with the above mass range, SNe la should 
have various lifetime depending on the main-sequence life- 
time of a secondary star. Therefore, SNe la are expected 
to continuously occur for a long time once they take place 
in the most massive binary system (e.g., Greggio & Ren- 
zini 1983). Moreover, most recently a new SNe progenitor 
scenario in which SNe la rate depends on the metallicity 
was suggested and it was shown that this scenario is suit- 
able for the explanation of the observed chemical evolution 
in the solar neighborhood (Kobayashi et al. 1998). How- 
ever, we fix the lifetime of SNe la to 1.5 Gyr for simplicity, 
following a result of Yoshii, Tsujimoto, & Nomoto (1996), 
who presented that the mean lifetime of the progenitors of 
SNe la is about 1.5 Gyr, based on the chemical evolution 
analysis in the solar neighborhood. 



2.2.3. The Energy Feedback 

One of the most difficult and most critical processes 
to model in galaxy formation simulations is the way in 
which the feedback from SNe affects the surrounding gas. 
Unfortunately, there is no well understanding of how it 
should be modeled. In a pioneering SPH simulation by 
Katz (1992) the energy feedback was implemented as de- 
position of purely thermal energy into the surrounding gas. 
Katz (1992) founds that this form of energy feedback is 
rather inefficient; the high densities of typical star- forming 
regions imply short cooling timescales that minimize the 
hydrodynamical effects of the feedback energy input. In 
real star-forming systems SNe are largely responsible for 
the multi-phase structure of the interstellar medium (e.g., 
Mckee & Ostrikcr 1977). Although they are associated 
with the cool and dense gas clouds, they can deposit much 
of their energy in a hot and low density wind. A sem- 
inal attempt at representing this multi-phase interstellar 
medium (ISM) has been carried out by Yepes et al. (1997) 
and recently Hultman & Pharasyn (1999) have adapted 
the Yepes et al. model to SPH simulations. However, they 
had to make a number of assumptions about the physics 
of the ISM and their models are not truly multi-phase, be- 
cause different phases are treated as a dynamically single 
entity. Another implementation assumed that the feed- 
back region around newly born stars evolves adiabatically 
until the multiple SNe II phase ends (e.g., Mori et al. 1997; 
Thacker & Couchman 2000). Although this model is 
somewhat artificial, it can provide strong effects on the 
formation of dwarf galaxies (Mori et al. 1997). 

While a number of models for energy feedback are pro- 
posed by previous authors as mentioned above, we adopt 
here the model proposed by Navarro & White (1993). 
This model assumes for simplicity that the energy pro- 
duced by SNe affects only the temperature and velocity 
field of the surrounding gas and its effect is implemented 
by increasing the internal and kinetic energy of the gas 
neighbors of each star particle by amount corresponding 
to the energy released by SNe. Thus, when i-th new born 
star particle ejects the energy of ^sn,« and deposits the 
energy of AEgNj m the J-th neighbor gas particle, the 
velocity of the j-th gas is altered by Avsn.j, which is cal- 



culated from 



fvEsN.i — Winj^hij) = m J (v j -Av SN . j + -Avg N:j ) 



(7) 

where Vj is the relative velocity to the i-th star, AugN j is 
a velocity perturbation directed radially from the i-th star, 
and /„ is an input parameter which controls the fraction 
of the available energy to perturb the gas velocity field. 
The rest of the energy of the SNe, (1 — f v )AEsw,j, con- 
tributes to increase the internal energy of j-th gas. When 
the term of Vj ■ Avsnj is assumed to be zero approximately 
(Navarro & White 1993), the magnitude of the Avsnj is 
then simply 

This model is consistent with the model of Katz (1992), 
in which f v — 0. The amount of effects of the energy feed- 
back on forming galaxies depends quite strongly on the 
value of /„. It is known that the kinetic feedback affects 
the history of star formation more strongly than the ther- 
mal feedback which quickly dissipates due to the radiative 
cooling in the high density region where stars can form 
(Navarro & White 1993). Thus, /„ controls the strength 
of the effect of SNe. 

We assume that each SN yields energy of e$N x 10 51 
ergs. According to high resolution ID simulations of SN 
remnant in Thornton et al. (1998), 90 % of an initial SN 
energy is lost in radiation in its early expansion phase, 
which is not resolved in our simulations. However, an ini- 
tial SN energy has not been established quantitatively yet. 
Therefore, we consider the available SN energy as a free 
parameter. 

In our code, there are two parameters to control the 
strength of the effect of SNe. To examine the effect of the 
difference in the strength of the energy feedback, we carry 
out the simulations adopting two extreme sets of these 
parameters, namely (e S N,fv) = (0.1,0) and (e S N,fv) = 
(4,0.9). The former means a minimum feedback model, 
corresponding to the result of Thornton et al. (1998) when 
an initial SN energy 10 51 ergs is assumed. The latter pro- 
vides a feedback strong enough to reproduce the observed 
CMR as shown later. 

3. THE MODEL 

Using the above code, we calculate the following semi- 
cosmological model which is almost the same as the one 
in K99. We consider an isolated sphere, as a seed galaxy, 
on which small-scale density fluctuations corresponding to 
a CDM power spectrum are superimposed. Here, we use 
Bcrtschinger's software COSMICS (Bertschinger 1995) in 
generating initial density fluctuations. To incorporate the 
effects of fluctuations with longer wavelengths, the density 
of the sphere has been enhanced and a rigid rotation cor- 
responding to a spin parameter, A, has been added. The 
initial condition of this model is determined by the fol- 
lowing four parameters: A, M tot , cs.in, and z c . The spin 
parameter is defined by 

7\E\ 1 / 2 

*=^67a. ( 9 ) 

where J is the total angular momentum of the system, E is 
the total energy, and M tot is the total mass of this sphere, 
which is composed of dark matter and gas; ag,in is the rms 
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mass fluctuation in a sphere of radius 8 h~ 1 Mpc, which 
normalizes the amplitude of the CDM power spectrum; z c 
is the expected collapse redshift. If the top- hat density 
perturbation has an amplitude of Si at the initial redshift, 
Zi, we obtain z c = 0. 36(5^(1 + Zi) — 1 approximately (e.g., 
Padmanabhan 1993). Thus, when z c is given, Si at Zi is 
determined. 

K99 found that the seed galaxy which has a slow rota- 
tion corresponding to A = 0.02 and the small-scale density 
fluctuations evolves into an elliptical-like system using the 
similar numerical simulation to that of this paper. This 
spin parameter is close to the minimum value possible in 
the CDM universe, according to the results of N-body sim- 
ulations (Barnes & Efstathiou 1987; Warren et al. 1992). 
In this paper we adopt this elliptical galaxy formation sce- 
nario. The main purpose of this paper is to study the evo- 
lution of seed galaxies with different masses. In addition 
we consider that the strength of feedback is crucial to give 
rise to mass-dependent evolution process. Thus, we simu- 
late the evolution of the seed galaxies with the total mass 
between 8 x 10 12 and 1 x 10 11 M Q , using the two codes 
with different strength of the SNe feedback as mentioned 
above. Parameters in each model are summarized in Ta- 
ble 2. Each model is labeled, for example, model 'Al', 
according to the total mass and the strength of the SNe 
feedback. Models A have (esN,fv) = (0-1,0), whereas 
(csn, fv) = (4, 0.9) in models B. The number in the model 
name represents the total mass of the seed galaxy. Here- 
after we refer to the mass sequence from model Al (Bl) to 
model A4 (B4) simply as "model A (B)" . In all the models, 
we fix as, in = 0.5 and z c — 3.5. 

Our simulations assume a flat universe (fi = 1) with 
a baryon fraction of fib = 0.1 and a Hubble constant of 
i?o = 50 km s _1 Mpc -1 . We carry out each simulation 
using 9171 particles for gas and dark matter respectively. 
The mass and spatial resolutions for each model are shown 
in Table 2. We simulate the evolution of each model from 
Zi = 40 to z = 0. The morphological evolutions of all the 
models are similar to the evolutions seen in Figure 1 of 
K99 and K01. Finally, nearly spherical stellar systems are 
formed at z = in all the models. 

4. DATA ANALYSIS 

To compare simulation results with the observational 
data of nearby elliptical galaxies directly, we have to de- 
duce photometric properties of simulation end-products at 
z = for various models. We use the same procedures as 
K01. Here we briefly describe the procedure to derive the 
photometric properties from simulation results. For this 
purpose, we show the analysis for model A2 as an example. 

In our simulations, the stellar particles contain the in- 
formation about their age and metallicity due to the self- 
consistent calculation of the chemical and dynamical evo- 
lution. By means of the population synthesis, we can de- 
rive the photometric properties of the stellar system from 
this information. Here, the spectral energy distribution 
(SED) of each stellar particle is assumed to be that of a 
single stellar population (SSP) that means a coeval and 
chemically homogeneous assembly of stars. Since the ob- 
servational data with which our results should be com- 
pared provide the luminosity distribution projected to a 
plane, we have to derive the projected distribution of SED 
from the three dimensional distribution of stellar parti- 



cles. Finally, we obtain the projected images as shown in 
Figure 5 of K01. Then the flux of each stellar particle 
is smoothed using a gaussian filter with the filter scale of 
1/4 of the softening length of the stellar particle. These 
images provide quite similar information to the imaging 
data obtained in actual observations. Thus we can obtain 
various photometric properties from these images in the 
same way as in the analysis of observational imaging data. 
In the following analysis, we use the images similar to the 
one displayed in Figure 5 of K01, but employing a 1001 x 
1001 pixel mesh to span the squared region with 100 kpc 
on a side. 



model A2 model B2 




0.1 1 10 0.1 1 10 

Radius (kpc) Radius (kpc) 

Fig. 1. — The B (upper panel) and V (lower panel) band surface 
brightness profiles in the x — y projections for models A2 (left) and 
B2 (right). The solid lines denote the r 1 /™ law best-fit for the data 
plotted as the solid symbols. The fitting parameters are shown in 
the lower left corner of each panel. The error bars are shown only 
for the solid symbols and correspond to the standard deviation (see 
text for details). 

In this paper, we use the data of SSPs of Kodama & 
Arimoto 97 model (Kodama 1997; Kodama & Arimoto 
1997). Kodama & Arimoto 97 model supplies the database 
of SSPs with two types of IMF: (x, M u , Mi) = (1.35, 60, 
0.1) and (x, M u , Mi) = (1.1, 60, 0.1) in the definition of 
Section 2.1. We adopt the data of SSPs with the IMF of 
(x, M u , M) = (1.35, 60, 0.1), while we use the IMF of (x, 
M u , Mi) = (1.35, 60, 0.4) in the numerical simulations. 
We confirmed that the photometric properties presented 
in this paper are hardly changed, when we employ another 
dataset of SSPs with the IMF of (x, M u , M) = (1.1, 60, 
0.1). A shallower IMF has a similar fraction of massive 
stars to the Salpeter IMF (x = 1.35) whose lower mass 
limit, Mi, is a little higher. Therefore, our use of the IMF 
of (x, M u , Mi) = (1.35, 60, 0.1) is justified. 

The left panels of Figure 1 shows the surface-brightness 
profiles for model A2 in the B and V bands obtained by 
setting the annuli of various radii in each band image. 
The error bars show the standard deviation in each annu- 
lus. We set the width of each annulus to half the soften- 
ing length. Since the error is calculated in the flux, like 
er 2 B =< fg > — < Jb > 2 , the error in the magnitude is 
defined, using the Taylor expansion, as 

° Mb = 2 ' 5 (<//> B ln(10) " 2<//>Mn(10)) ' (10) 
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These surface brightness profiles are in excellent agreement 
with the Sersic (r 1 /™) law (Sersic 1968), 

J(r) = I e lQ{- b ^/r^-l]}^ (n) 

where we adopt b n = 0.868n — 0.142, so that the effec- 
tive radius, r e , equals the half light radius in the range 
0.5 < n < 16.5 (Caon, Capaccioli, & D'Onofrio 1993) and 
I e is the surface brightness at r e . This corresponds to the 
de Vaucouleurs (r 1 ' 4 ) law (de Vaucouleurs 1948), when 
n = 4. The solid lines in Figure 1 show the best-fits in 
applying equation (11). In fitting, we use the data with 
Hb < 27 mag arcsec~ 2 irrespective of the observed band. 
The best-fit parameters are shown in the lower left corner 
of each panel in Figure 1. 



model A2 model B2 

LO 





-0.5 0.5 1 1.5 -0.5 0.5 1 1.5 

Log(r) Log(r) 

Fig. 2. — Color (top), luminosity weighted metallicity (2nd), lu- 
minosity weighted age (3rd), and [Mg/Fc] (bottom) gradients in the 
x — y projection for models A2 (left) and B2 (right). The solid lines 
show the best-fit linear relations for the data plotted as the solid 
symbols. The gradients in best-fit lines are shown in the lower left 
corner of each panel. The error bars are shown only for the solid 
symbols in top and 2nd panels and for the data of r < 5 kpc in 3rd 
and bottom panels. 

The left panels of Figure 2 shows the color, metallicity, 
and age gradients for model A2. The color gradients are 
obtained by setting annuli of various radii in the B and 
R band images and subtracting the R band magnitude 
from the B band magnitude in each annulus. We made 
the projected images also for the metallicity and age, and 
obtained their radial profiles. Points indicate the average 
value in each annulus in the x-y projections. The width 
of each annulus is set to half the softening length. The 
error bars show the standard deviation in e ach annulus . 
The error in the color is written as <jb— r — \J &m b T ' 

The profiles of color and metallicity [defined as \og(Z / Zq)] 
are fitted by linear regression. In fitting, we excluded the 
data at radii less than the softening length and greater 
than the radius at which the B band surface brightness is 
[iB = 24.5 mag arcsec -2 (open symbols in Fig. 2), because 



the inner region is affected by the smoothing of gravita- 
tional forces whereas the number of particles within an 
annulus is too small in the outer region. The best-fit gra- 
dients are shown in the lower left corner of each panel in 
Figure 2 and also in Table 3. The data obtained from the 
simulation results have large errors leading to large errors 
in the calculated gradients. However negative gradients 
are clearly seen in colors and metallicity, i.e., the color 
(metallicity) at the center is redder (higher) than that 
in the outer region. Although we show only a luminos- 
ity weighted metallicity to compare with the observational 
data, the metallicity gradient weighted by the luminosity 
is slightly steeper than that weighted by the mass in all 
the models. It is known that a typical elliptical galaxy has 
A(B - i?)/Alog(r) = -0.09 ± 0.02 (Peletier et al. 1990) 
and Alog(Z/Z )/Alog(r) = -0.30 ±0.12 (Kobayashi & 
Arimoto 2000, and references therein). Although the gra- 
dients obtained from the simulation results have large er- 
rors, the median gradients of both color and metallicity 
for model A2 are steeper than those observed. We can 
see that the color gradients are caused by the metallicity 
gradients, because the age is equally old irrespective of the 
radius. The radial profile of the age has a slightly positive 
gradient. This positive gradient means that star forma- 
tion persists over longer time in the inner region than in 
the outer region, due to dissipative infall of gas. Since the 
residual gas is polluted by past star formation, metal rich 
stars are formed in the central region where the duration 
of star formation is long (K99). This is the origin of the 
color and metallicity gradients. Recent observations using 
the metal absorption lines provide no (for ellipticals) or 
positive (for SOs) gradients of the age (Kuntschner 1998). 
Since the age gradient in Figure 2 is too shallow, this gra- 
dient seems not to be detected in the current observations 
using the absorption lines. 




0.5 1 2 5 10 -24 -22 -20 -18 



R h ,v 0/n) M fcv (1/n) 

Fig. 3. — Comparison of the effective radius (left) and total 
magnitude (right) derived by different methods. The upper panels 
compare the value derived from the r 1 /™ law fitting in the x-y pro- 
jection with that in the x-z projection. The lower panels compare 
the value derived from the r 1 /™ law fitting (square) with the value 
derived from the large aperture photometry (circles) in the x-y pro- 
jection (see text for details). The open (solid) symbols indicate the 
value for model A (B). 

We show the surface-brightness profiles and the color, 
metallicity, and age gradients for model B2 in the right 
panels of Figures 1 and 2 in order to evaluate the effect of 
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the difference in the strength of the SNe feedback on these 
profiles and gradients. We analyzed these profiles and gra- 
dients also for the other models. Table 3 summarizes the 
results. In Table 3, we can see that the strong feedback of 
model B leads to a small n in surface-brightness profiles 
and a shallow slope in color and metallicity gradients, com- 
pared to the ones of model A. Moreover, with decreasing 
mass of the system, n decreases and the color and metal- 
licity gradients becomes shallower in model B. In other 
words, the strong feedback affects these profiles and gra- 
dients more severely in lower mass systems. Similar results 
are obtained in Carlberg (1984a); Mori, Yoshii, & Nomoto 
(1999). The mass dependence of n is consistent with the 
correlation between n and luminosity in the observed ellip- 
tical galaxies qualitatively (Caon, Capaccioli, & D'Onofrio 
1993). On the other hand, there is no correlation between 
the observed color and metallicity gradients and the lumi- 
nosity in luminous elliptical galaxies (Peletier et al. 1990; 
Kobayashi & Arimoto 2000, and reference therein), and 
in less luminous elliptical galaxies (Mb > —19) there is 
still a large uncertainty (e.g., Vader et al. 1988). 

The global photometric properties, such as the total lu- 
minosity, the colors, and the effective radius, are also ob- 
tained from the projected image data. In this paper, we fo- 
cus on the CMR, the Kormendy relation, and the [Mg/Fe] - 
magnitude relation. Thus we have to analyze the total 
magnitude, the color, the effective radius, and [Mg/Fe] us- 
ing the projected images. The total magnitude and the 
effective radius can be derived in two different ways. One 
method is to fit the surface brightness profile to the r 1 /" 
law (eq. [11]). This fitting formula directly provides the 
effective radius, and the integration from the center to the 
effective radius gives half the total luminosity. The other 
method is to assume that the aperture photometry with 
a sufficiently large aperture provides the total value. We 
assume that the luminosity derived with the 99 kpc aper- 
ture is the total luminosity. Then the effective radius is 
defined as the radius within which half the total luminos- 
ity is contained, and calculated by counting the luminosity 
from the center in the projected images. Throughout this 
paper we set the center to the position of a pixel which 
has the maximum V band luminosity. Figure 3 compares 
the effective radius and the total magnitude derived by the 
different methods. These panels show that there is no sys- 
tematic difference. Moreover Figure 3 indicates that the 
total magnitude and the effective radius do not depend on 
the direction of the projection. It means that the end- 
products are nearly spherical and the luminosity profiles 
are well approximated by the r 1 /™ law (see also Figs. 1). 
Finally, we use the former method in the following discus- 
sion. Table 3 shows the total magnitude in the B, V , and 
K bands and the effective radius in the V and K bands 
for each model. These data are derived from the fitting of 
the luminosity profiles in the x-y projections and used in 
the next section. 

In observation of elliptical galaxies, colors and [Mg/Fe] 
are usually obtained using a fixed aperture. Table 4 shows 
the U—V and V — K colors and [Mg/Fe] for each model. In 
the next section, we will compare the colors and [Mg/Fe] 
of the simulation end-product in each model with the 
observed values measured by BLE92a and Trager et al. 
(2000a) . The values in Table 4 have been obtained by the 



analysis using the same aperture as BLE92a and Trager 
et al. (2000a). We made the projected image also for the 
luminosity weighted abundance of Mg and Fe, and ana- 
lyzed the ratio of the mean abundance of Mg to that of Fe 
within the aperture. 
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Fig. 4. — The comparison of V/cr for the simulation end-products 
and for galaxies in the Virgo and Coma clusters (small dots, Davies 
et al. 1983). The triangles (circles) connected by dotted (solid) 
lines indicate the V/cr of model A (B). 

Due to calculating dynamical evolution, we can analyze 
the dynamical properties of end-products. K01 showed 
that the elliptical galaxy formation model used in this pa- 
per leads to a small ratio of the rotational velocity to the 
velocity dispersion (V/ a) , which is consistent with the ob- 
served V/a in bright elliptical galaxies. The observed V/a 
in elliptical galaxies increases with decreasing luminosity 
(Davies et al. 1983). Now we consider the galaxies with 
different masses, we analyzed V/a for all the models, to 
examine the mass dependence of V/a. Figure 4 shows 
V/a against the B band absolute magnitude. In model A, 
V/a increases slightly with decreasing luminosity, like the 
tendency of the observational data. On the other hand, 
the strong feedback of model B seems to lead to a low 
V/a, compared to the one of model A, in low mass sys- 
tems. However, V/a of simulation end-products could not 
be measured accurately, because their rotation curves do 
not provide a clear flat shape, owing to their slow rota- 
tion and a poor spatial resolution of numerical simulation. 
In other words, V/a indicated in Figure 4 has ambigu- 
ity. Thus we conclude that all the models provides a small 
V/a around 0.1 ~ 0.5 irrespective of the luminosity or the 
feedback model. Our model cannot explain the existence 
of less luminous elliptical galaxies with a large V/a, es- 
pecially in model B, tough the observed (V/a)s for less 
luminous galaxies have a broad range. We suppose that 
these systems might be formed from another initial con- 
dition, e.g., with a larger A or a smaller <78,i n (K01). This 
analysis made it clear that the final stellar systems which 
we focus on have a small V/a and are supported mainly 
by random motions in terms of kinetics. In the following 
sections, we concentrate on the photometric properties. 

5. GLOBAL SCALING RELATION 

The photometric properties of the simulation end- 
products at z = which are derived by the data analy- 
sis described in Section 4 are compared with the observed 



8 



D. Kawata 



global scaling relations of elliptical galaxies in this section. 
We first compare the simulation results with the observa- 
tional data in the color versus magnitude diagram. Next 
we discuss the sizes of the simulation end-products, com- 
paring them with the observed Kormcndy relation. Finally 
we examine the correlation between the total magnitude 
and the abundance ratio of Mg to Fe, i.e., the [Mg/Fe]- 
magnitude relation, which is sensitive to the star formation 
history. 
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(BLE92b). This gives the luminosity distance of 87.1 Mpc 
for the Coma. We assume that the angular diameter dis- 
tance equals the luminosity distance, because the redshift 
of the Coma cluster (z ~ 0.023) is nearly zero cosmologi- 
cally. Then the aperture size of 11 arcsec at the distance 
of the Coma cluster corresponds to ~ 5 kpc. The data for 
simulation end-products are tabulated in Tables 3 and 4. 
Although model A has a tendency that smaller galaxies 
have a bluer color, the slope is too shallow to reproduce 
the observed slope of the CMR in the Coma cluster. On 
the other hand, the slope of model B is steeper than that of 
model A, and model B well reproduces not only the slope 
but also the zero-point of the observed CMR. The suc- 
cess of model B is owing to the action of mass-dependent 
galactic wind as shown later in detail. 
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Fig. 5. — The comparison of the CMRs for the simulation end- 
products and the Coma cluster galaxies (small dots, BLE92a) in 
the aperture of 5 kpc. The triangles (circles) connected by dotted 
(solid) lines indicate the CMRs of model A (B). The dashed line 
shows the CMR fitted to the Coma cluster galaxies (BLE92b). The 
cross denotes the results of the same model as model A2 but with 
c* = 1 (see Section 2.1). 



5.1. The Color-Magnitude Relation 

We construct the CMRs in model A and B. Comparing 
model A with model B, we study the effect of the differ- 
ence in the strength of the SNe feedback. Figure 5 shows 
the comparison of the simulation end-products and the 
observed galaxies in the Coma cluster in the V — K and 
U — V CMR. The data for galaxies in the Coma cluster 
are the observed CMR of BLE92a. Since there is no dif- 
ference between SOs and ellipticals in the scaling relations 
which we discuss, we do not distinguish SOs from ellipti- 
cals. BLE92a supplies the U — V and V — K colors which 
refer to an aperture size of 11 arcsec and the V band to- 
tal magnitude derived from a combination of their data 
and the literature. Throughout this paper we adopt the 
distance modulus of the Coma cluster of m — M = 34.7 
mag; the Virgo distance modulus is m — M — 31.01 (Gra- 
ham et al. 1999) and the relative distance modulus of 
the Coma with respect to the Virgo is m — M = 3.69 
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Fig. 6. — The CMRs for the simulation end-products in the 99 
kpc aperture. The triangles (circles) connected by dotted (solid) 
lines indicate the results of model A (B). The dashed lines are the 
same as that in Figure 5. 

The slope of model A is just an aperture effect. El- 
liptical galaxies have a color gradient that the color at 
the center is redder than that in the outer regions (e.g., 
Peletier et al. 1990) and the end-products in our simula- 
tions also have gradients (Fig. 2 in Section 4). Thus the 
color within a fixed aperture provides the color in a more 
central region for larger galaxies, and the color of large 
galaxies becomes redder than that of small galaxies even 
if the mean color of the whole galaxy is the same between 
the large and small galaxies. Figure 6 shows the CMR 
for the mean color within the aperture of 99 kpc, which 
covers almost the whole galaxy in all the models which we 
examined. The mean color of the whole galaxy in model 
A stays nearly constant when the mass varies. Hence the 
slope of model A in Figure 5 is caused only by the small- 
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ness of the aperture size. On the other hand, model B 
shows a significant slope even when a large aperture size 
is specified as shown in Figure 6. Aperture effect cannot 
be ignored when we discuss the CMR observed in an aper- 
ture of a small size. Unfortunately, accurate colors eval- 
uated within various apertures have never been provided 
from any observations. Recently Scodeggio (2001) showed 
that the CMR derived using color measurements within 
the effective radius is significantly flatter than those based 
on fixed-aperture color measurements. Unfortunately, the 
photometric accuracy of their data was not good enough 
to derive an accurate CMR. We would like to stress that 
an observation in apertures of various sizes is quite im- 
portant in understanding not only the CMR but also the 
color gradients in elliptical galaxies. 
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Fig. 7. — The mctallicitics (upper panel) and ages (lower panel) 
against the absolute V band magnitude for each model. The trian- 
gles (circles) connected by dotted (solid) lines indicate the results of 
model A (B). The open (solid) symbols denote the values evaluated 
in the 5 kpc (99 kpc, spread over almost the whole galaxy) aperture. 

There are two equally plausible interpretations of the 
CMR (Worthcy 1994). The sequences of colors and line 
strengths among elliptical galaxies can almost equally well 
be attributed to either age difference or metallicity differ- 
ence. In the simulation this degeneracy can be broken 
completely, because we can directly analyze the metallic- 
ity and age in the simulation end-products. Figure 7 shows 
the metallicity and age for each model. Here the metal- 
licity and age mean the luminosity weighted value. The 
age is very old irrespective of the models and the aper- 
ture. The metallicity slope exhibits a similar behavior to 
the color as a function of the luminosity. Model A has a 
weak slope in the metallicity-magnitudc diagram of Fig- 
ure 7 when the aperture size is 5 kpc, though the slope 



almost vanishes when the 99 kpc aperture is used. On the 
other hand, model B shows a significant slope irrespective 
of the aperture size, although a large aperture size leads 
to a shallower slope. Thus it is concluded that the slope 
of CMR in Figure 5 is caused solely by the effect of the 
metallicity. 
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Fig. 8. — Time variations of the event rate of SNe II (normal 
lines) and SNe la (thick lines) for all the models. 

As a result, model B succeeds in reproducing the CMR 
by causing a decrease in the metallicity of the stellar con- 
tents for low-mass systems. On the other hand, model A 
completely fails. To see the origin of this difference, we 
examine the star formation history for each model. Fig- 
ure 8 shows the histories of SNe II and SNe la event rates 
normalized by the total mass of the system. In our sim- 
ulations, the histories of SNe II events completely trace 
those of star formation, because we assume the instanta- 
neous recycling (see Section 2.1 for details). On the other 
hand, SNe la occur 1.5 Gyr behind SNe II. We can clearly 
see that in model B star formation abruptly ceases around 
t = 2.5 ~ 3 Gyr (z = 2 ~ 1.7). The dynamical evolu- 
tion of the system clearly demonstrates that this is caused 
by the so-called galactic wind. Figure 9 which shows the 
morphological evolution of model B4. Around z ~ 1.7, gas 
particles begin to blow out from the stellar system, i.e., the 
galactic wind occurs. All the gas particles overcome the 
binding energy of the dark matter and stars, and escape 
from the system. In addition we found that the mass frac- 
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Fig. 9. — Time evolution of the system in model B4. The upper and lower panels shows the distributions of stellar particles and gas 
particles respectively. Upper (lower) panels show the x-z projection of the particles and measure 15 (150) kpc across. 



tion of the ejected gas increases with decreasing mass of 
the system in model B, as seen in Figure 10. The ejected 
gas mass in Figure 10 is defined as the total mass of all the 
gas particles whose galactocentoric radius is greater than 
20 kpc. Since the ejected gas cannot contribute to the fur- 
ther metal enrichment, a higher mass system suffers from 
more enhancement in heavy elements. This mass depen- 
dence of the mass fraction of the gas ejected by the galac- 
tic wind is the origin of the metallicity-magnitude relation 
shown in Figure 7, which in turn causes the CMR shown 
in Figure 5. There is little difference models Bl and B2 in 
Figure 10. In this luminosity range, the aperture effect is 
dominant factor in explaining the observe slope. However, 
the metallicity-magnitude relation caused by the galactic 
wind is essential in order to reproduce the global CMR in 
the observed luminosity range. 

It is notable that the galactic wind occurs immediately 
after SNe la ignite. It is clearly shown that SNe la are 
important for the galactic wind. The event rate of SNe la 
exceeds that of SNe II after the peak of star formation. 
SNe la occur when the gas is exhausted by star forma- 
tion and the gas density becomes low. Thus SNe la can 
strongly affect the evolution of the residual gas. 

On the other hand, Figure 10 shows that almost all the 
gas is exhausted by z — irrespective of the mass of the 
system in model A. As a consequence, model A fails to 
reproduce the observed CMR. The difference in results 
of models A and B is the strength of the effect of SNe. 
As a result, a strong effect of the SNe feedback like the 
one included in model B is required to explain the ob- 



served CMR. Moreover we can clearly see from the results 
of model B that the mass fraction of the gas ejected by 
the galactic wind depends on the mass of the system and 
that SNe la play a quite important role in the outbreak of 
galactic wind. 
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Fig. 10. — The ratio of the ejected gas mass (M g qw) at 2 = 
to the initial gas mass (M g ; n ;). The triangles (circles) connected 
by dotted (solid) lines indicate the data of model A (B). Here, the 
ejected gas is defined as the sum of all the gas particles whose galac- 
tocentoric radius is greater than 20 kpc. 



5.2. The Kormendy Relation 



The Chemo-Dynamical Evolution of Elliptical Galaxies 



11 



Numerical simulations have a strong advantage over 
pure chemical evolution studies and semi-analytic models 
of galaxy formation in that simulations can provide infor- 
mation about the structure of the end-products. Owing to 
this advantage, we can discuss the scaling relation which 
prescribes the size of the elliptical galaxies, i.e., the Kor- 
mendy relation. Figure 11 displays the comparison of the 
Kormendy relations for the simulation end-products and 
the Coma cluster galaxies both in the V and K bands. 
We refer to the data of the Coma cluster galaxies of Pahre 
(1999). When we derive the absolute magnitude and the 
effective radius in the kpc unit from the data set in Pahre 
(1999), we assume the same distance modulus as men- 
tioned above. In Figure 11, both models have a tendency 
that higher mass galaxies have larger effective radii, which 
is qualitatively consistent with the Kormendy relation for 
the observed elliptical galaxies. However the slope is too 
shallow in both models to consistent with the observed re- 
lation quantitatively. In particular, the slope of model B 
is much shallower than that of model A. This is due to the 
expansion of the low mass systems induced by the strong 
feedback. It is concluded that the strong feedback which 
is required to explain the observed CMR is unsuitable for 
explaining the Kormendy relation. The failure in repro- 
ducing the Kormendy relation poses a serious problem to 
the galactic wind scenario. 
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Fig. 11. — The comparison of the Kormendy relations for the 
simulation end-products and the Coma cluster galaxies (small dots, 
Pahre 1999) in the V (upper panel) and K (lower panel) bands. 
The triangles (circles) connected by dotted (solid) lines indicate the 
data for model A (B). The cross denotes the results of the same 
model as model A2 but with c = 1 (see Section 2.1). 



5.3. The [Mg/Fe]-Magnitude relation 



Our numerical code takes account of both SNe II and 
SNe la, and follows the evolution of the abundances of sev- 
eral chemical elements (C, O, Ne, Mg, Si, and Fe). There- 
fore we can examine the abundance ratio of chemical ele- 
ments in the simulation end-products. Observational data 
indicate that [Mg/Fe] correlates with the luminosity, i.e., 
the [Mg/Fe]-magnitude relation (see also Section 1). We 
compare the [Mg/Fe]-magnitude relation for the simula- 
tion end-products with the observational data of local field 
and group elliptical galaxies. We use the observational 
data of [Mg/Fe] in Trager et al. (2000a) and the abso- 
lute B band magnitudes in Trager et al. (2000b). Trager 
et al. (2000a) derived [Mg/Fe] in the r e . B /2 and r e>B /8 
apertures, using an extension of the Worthey (1994) mod- 
els that incorporates nonsolar line-strength response func- 
tions by Tripicco & Bell (1995). We use the data for 
model 4 of Trager et al. (2000a), which is the most suc- 
cessful model according to their paper. In Figure 12 we 
compared the r e .s/2 aperture data with [Mg/Fe] for the 
simulation end-products in the r e ,s/2 aperture. Since all 
the models do not reproduce the observed effective radius 
as shown in Section 5.2, the physical aperture size in all 
the models is not consistent with that applied in the ob- 
served galaxies with the same luminosity. Fortunately, the 
observed gradient of [Mg/Fe] is quite weak in these sam- 
ples according to Trager et al. (2000a). On the other 
hand, Figure 2 shows the [Mg/Fe] gradients for models A2 
and B2. Although there is a little gradient of [Mg/Fe] in 
the simulation results of model A, the gradient becomes 
shallower within the effective radius irrespective of models. 
Thus the difference in the aperture size is not so impor- 
tant in this case. The observed data show a clear tendency 
that [Mg/Fe] increases with the galactic luminosity. How- 
ever both models A and B are incapable of reproducing 
this tendency. In Figure 8 the star formation histories in 
both models A and B are nearly identical for the whole 
mass range of the system, although there is a clear differ- 
ence between models A and B. Consequently, [Mg/Fe] is 
constant irrespective of the luminosity of the system. In 
addition [Mg/Fe] in model B is systematically higher than 
that in model A. The reason is that in all of model B the 
galactic wind leads to the cessation of star formation be- 
fore the metal enrichment by SNe la progresses, while in 
all of model A star formation continues after SNe la ignite. 

Model B provides a little different result in [Mg/Fe]- 
magnitude diagram from what pure chemical evolution 
studies of the galactic wind scenario predicted (e.g., Mat- 
teucci 1994). In pure chemical evolution studies, the ef- 
ficiency of star formation is assumed to be constant irre- 
spective of the mass of the systems or an increasing func- 
tion with decreasing mass. Therefore, the galactic wind in 
higher mass systems need to occur later than that in lower 
mass systems in order to explain the observed CMR, and 
the duration of star formation is longer in the higher mass 
systems than that in the lower ones. The elements of Mg 
and Fe are mostly produced by SNe II and SNe la respec- 
tively, and SNe la have a longer delay than SNe II with 
respect to the formation of stars. A galaxy with a longer 
time duration of star formation is much more enriched 
by SNe la and gets a lower [Mg/Fe]. Hence pure chemi- 
cal evolution studies predict that [Mg/Fe] is a decreasing 
function of the galactic mass and luminosity, which pro- 
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vides an opposite slope to the observation (Matteucci & 
Tornambe 1987). On the other hand, in model B [Mg/Fe] 
is constant irrespective of the luminosity. Comparison of 
the peak SNe II event rate for each model in Figure 8 
shows that the efficiency of star formation increases with 
mass of the system. It means that the strong feedback 
causes a self-regulation of star formation in lower mass 
systems (see also Lia, Carraro, & Salucci 2000). Con- 
sequently, the galactic wind in all the models occurs at 
nearly the same time, though the higher mass systems con- 
verts a larger fraction of gas into stars as seen in Figure 
10. Therefore, our numerical simulation which calculates 
the dynamical evolution self-consistently reveals that the 
mass-dependent self-regulation of star formation is caused 
by the strong feedback, so that the problem of the galactic 
wind scenario in the [Mg/Fe]-magnitude relation is im- 
proved. Finally, the high constant [Mg/Fe] in model B 
is caused by both effects of the cessation of star forma- 
tion in an early epoch irrespective of the mass of systems 
and the mass-dependent self-regulation of star formation, 
which are induced by the feedback strong enough to re- 
produce the observed CMR. However, the observed slope 
is still not reproduced completely. 




Fig. 12. — The comparison of the simulation end-products with 
the early-type galaxies in the observed [Mg/Fe] (Trager et al. 2000a) 
vs. Mb diagram. The triangles (open circles) connected by dotted 
(solid) lines indicate the data for model A (B). The cross denotes 
the results of the same model as model A2 but with c* = 1 (see 
Section 2.1). 



6. DISCUSSION AND CONCLUSIONS 

We have numerically studied the dynamical and chem- 
ical processes of the formation of elliptical galaxies in 
the CDM universe, in order to examine the origin of the 
mass-dependence of the photometric properties of ellipti- 
cal galaxies. Our numerical simulations take account of 
both SNe II and SNe la, and follow the time evolution of 
the abundances of several chemical elements (C, O, Ne, 
Mg, Si, and Fe). In this paper we paid a special attention 
to the CMR of elliptical galaxies. In addition we examined 
the Kormendy relation, which prescribes the size of ellip- 
tical galaxies, and the [Mg/Fe]-magnitude relation, which 
provides a strong constraint on the star formation history. 
Because the SNe feedback is likely to be a crucial determi- 
nant of the properties of the end-products, we compared 
two different strengths of the SNe feedback, the minimum 



(model A) and strong (model B) feedback. Although our 
feedback implementations are very crude, we could exam- 
ine whether or not the SNe feedback should strongly affect 
the evolution of the system to explain the global scaling 
relations of elliptical galaxies. 

In conclusion, the strong influence of SNe like the one 
adopted in model B is required to reproduce the observed 
CMR as shown in Figure 5. We found that the feedback af- 
fects the evolution of lower mass systems more strongly, so 
that a larger fraction of gas is blown out in a lower mass 
system (Figure 10). Consequently, higher mass systems 
become more metal rich and have redder colors than lower 
mass systems (Figure 7) . These results are consistent with 
the predictions of pure chemical evolution studies based on 
the galactic wind scenario (e.g., Arimoto & Yoshii 1987). 
The mass dependence of the metal enrichment has already 
been examined qualitatively also in the realistic numeri- 
cal studies. However, previous authors adopted a simple 
model of elliptical galaxy formation in which an elliptical 
galaxy is formed by the monolithic collapse of a single ho- 
mogeneous proto-galactic cloud (e.g., Carlberg 1984a, b; 
Friaca & Terlevich 1998; Mori, Yoshii, & Nomoto 1999) 
or the mergers of two disk galaxies (e.g., Bekki & Shioya 
1998, 1999). We adopted the scenario that the ellipti- 
cal galaxies are formed through the clustering of the sub- 
clumps caused by the initial small-scale density perturba- 
tions. This modeling is based on the CDM cosmology and 
has a reasonable cosmological ground. Thus the present 
model is more realistic than those used by previous au- 
thors. We showed that the mass-dependence of metal en- 
richment is naturally induced in this scenario. The present 
work demonstrates that many conclusions of pure chemical 
evolution studies based on the monolithic collapse picture 
are valid also in elliptical galaxy formation through the 
hierarchical clustering. 

The assumed energy of each SN (4 x 10 51 ergs) in model 
B might seem to be too strong. It might suggest that an- 
other mechanism ignored in this study is needed in order 
to enhance the effect of the SNe feedback on the evolution 
of systems. For example, the UV background radiation 
is a plausible candidate, because it is possible to suppress 
the cooling in low mass systems, as discussed later. 

We were able to evaluate relative importance of two 
types of SNe, SNe II and SNe la. SNe II occur imme- 
diately after star formation events, while there is a delay 
of 1.5 Gyr between star formation and SNe la. Therefore, 
SNe II occur in dense gaseous environments still actively 
forming stars, whereas SNe la act on tenuous gas left after 
star formation. This difference makes SNe la the major 
trigger of the galactic wind. This result is a novel one ob- 
tained from our numerical simulations which calculate gas 
dynamics, radiative cooling, star formation, and the SNe 
II and SNe la feedback self-consistently, and means that 
the energy feedback of SNe la is of fundamental impor- 
tance for the evolution of elliptical galaxies. However, it 
should be noted that our SNe modeling is very simple. As 
mentioned in Section 2.2.2, we neglect the spread of the 
lifetime of SNe la and assume that all SNe la explode 1.5 
Gyr after star formation. This assumption might cause 
too great impact on the interstellar medium. Modeling of 
SNe II is also simplified by an assumption of the instan- 
taneous recycling. Our average time step was ~ 1.7 x 10 5 
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yr. It is considerably shorter than the typical lifetime of 
SNe II progenitors. The effects of feedback of SNe II might 
also be overestimated. In order to compare the effects of 
SNe II and SNe la quantitatively, it is required to adopt 
more sophisticated models of SNe la and SNe II which con- 
sider the realistic spread of lifetimes of SNe la and SNe II 
progenitors. We will address this issue in a forthcoming 
paper. Nevertheless, our results are still important in that 
it is clearly shown that the delayed SNe feedback, such 
as SNe la, affects the dynamical evolution of interstellar 
medium more efficiently than the instantaneous one. 

In Table 5 we compare the SNe II and la rates of the 
simulation end-products at z = to the observed rates 
at nearby galaxies (Caperrallo, Evans, & Turatto 1999). 
The present SNe rates are obtained by the average of the 
number of SNe which occurred within 1.5 Gyr and 2 Gyr. 
These values slightly depend on the duration of the sam- 
pling in model A, as expected from Figure 8. In all the 
models of model A SNe II rate is too high to reproduce the 
observed rate, which suggests few or no SNe II. Although 
model B is in agreement with no SNe II, the SNe la rate 
is inconsistent with the observation, which suggests a non- 
zero event rate. We anticipate that SNe la model which 
takes into account the spread lifetime of SNe la progeni- 
tors would eliminate disagreement between these SNe la 
rates, because the continuous SNe la is expected after star 
formation ceased. 

In analyzing the colors of simulation end-products, we 
payed attention to the aperture effect by taking the same 
aperture size as adopted in the actual observations. The 
color evaluated in an aperture fixed to a small size leads 
to a non-vanishing slope in the color-magnitude diagram, 
even if the sample galaxies have an identical mean color 
(i.e., the color averaged over the entire galaxy) as seen in 
model A. This is because elliptical galaxies generally have 
color gradients such that the center of the galaxy is redder 
than the outer regions. Hence we would like to stress that 
the aperture effect should not be ignored when we discuss 
the CMR observed in an aperture fixed to a small size. 

Our numerical models failed to quantitatively reproduce 
the Kormendy relation, which prescribes the size of ellip- 
tical galaxies. Both models A and B give a significantly 
shallower slope than the observed one. Especially, model 
B, which can reproduce the observed CMR, causes too 
large effective radii in low mass systems due to the ex- 
pansion of the system induced by the strong feedback. It 
means that the feedback strong enough to explain the ob- 
served CMR leads to not only the galactic wind but also 
the expansion of the system. This result rises a serious 
problem to the galactic wind scenario. 

The diagram of the abundance ratio [Mg/Fe] versus the 
galactic luminosity provides an important diagnostic to in- 
fer the star formation history in galaxies, due to the differ- 
ent timescales for the Mg and Fe production. The observed 
[Mg/Fe]— magnitude relation indicates that more luminous 
galaxies have larger [Mg/Fe]. However, both models A 
and B leads to nearly constant [Mg/Fe] as a function of 
the mass, and [Mg/Fe] of model B is substantially larger 
than that of model A due to the short duration of star 
formation. Both of models thus failed to reproduce the 
observed [Mg/Fe]-magnitude relation. However, we must 
bear in mind that theoretical prediction of [Mg/Fe] still 



remains ambiguous, because the ratio depends on some 
unknown physical processes, such as the lifetime of the 
progenitor of SNe la and the nucleosynthesis yields of SNe 
II and SNe la. 

This discrepancy in the [Mg/Fe]-magnitude relation has 
already been known as a problem of the galactic wind sce- 
nario. As mentioned in Section 5.3, our numerical sim- 
ulation in the strong feedback model, i.e., model B, re- 
duces the discrepancy, compared to the prediction of pure 
chemical evolution studies, due to the mass-dependent self- 
regulation of star formation. However, to explain the ob- 
served trend quantitatively, we still have to consider what 
physical process ignored in the present study can lead to a 
small [Mg/Fe] in lower mass systems. Faber, Worthey, 
& Gonzalez (1992) proposed three alternative scenar- 
ios to solve this discrepancy: i) different star formation 
timescales; ii) a selective loss of Mg and Fe; and iii) a se- 
lective production of Mg vs. Fe, i.e., variable IMF. There is 
no physical reason yet why the IMF depends on the mass 
of the system. It is a natural consideration that the IMF 
depends on the local physical condition of the interstellar 
gas rather than the mass of the host system. Thus the 
scenario iii) sounds infeasible. In our numerical simula- 
tions, ii) did not occur. As a physical process to realize 
i) we can consider the UV background radiation which we 
ignored in this paper. The UV background radiation can 
suppress the cooling and star formation more efficiently 
in lower mass galaxies (Efstathiou 1992), which is also 
confirmed by realistic numerical simulations (e.g., Quinn, 
Katz, & Efstathiou 1996; Thoul & Weinberg 1996; Wein- 
berg, Hernquist, & Katz 1997). Thus it is possible to 
decrease the efficiency of star formation and prolong star 
formation with decreasing mass of the system. A forth- 
coming paper will discuss the effect of UV background 
radiation on the [Mg/Fe] in lower mass galaxies. 

The galaxy formation model used here parameterizes 
the properties of the seed galaxy by the total mass, the 
spin parameter, the amplitude of the small scale density 
perturbations, and the collapse redshift, i.e., M tot , A, as,i n , 
and z c . This approach turned out to be quite useful in 
studying the evolution of a single seed galaxy in detail. 
In this paper we studied how the properties of the end- 
product depend on the total mass of the system, M tot , 
when other parameters are fixed. However, little is known 
about whether and how the above parameters depend on 
each other. On the other hand, in cosmological simula- 
tions, the above parameters which are assumed to be in- 
dependent of each other in this paper are automatically 
determined in the process of structure formation in the 
universe. Hence, a cosmological simulation which is able 
to spatially resolve the process of the formation and evo- 
lution of individual galaxies is indispensable in studying 
the connection between the evolution of galaxies and that 
of the universe. In future work we intend to investigate 
what is the origin of the observed scaling relations of the 
elliptical galaxies, i.e., what controls those relations, using 
a high-resolution cosmological simulation. 
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Table 1 

Nucleosynthesis products of SNe II and SNe Ia 



Element 


Synthesized 
Type II (M iiSN eii) 


mass a (M Q ) 

Type Ia (M ljSNcIa ) 


Solar abundance 13 


12 C 


7.91 x 10~ 4 


7.32 x 10" 5 


0.30 x 10~ 2 


16Q 


2.01 x lO" 2 


2.17 x 10~ 4 


0.96 x 10~ 2 


20 Ne 


2.24 x 10" 3 


3.06 x 10~ 6 


0.16 x 10~ 2 


24 Mg 


9.38 x 10~ 4 


1.29 x 10" 5 


0.52 x 10" 3 


28 Si 


1.10 x 10" 3 


2.27 x 10~ 4 


0.65 x 10" 3 


56p e 


7.79 x 10~ 4 


9.29 x 10~ 4 


0.13 x 10~ 2 



a The Salpeter IMF with M u = 60 M and Mi = 0.4 M, 
b Ferrini et al. 1992 



Table 2 
Model Parameters 



Model 




M tot 


Particle Mass (M Q ) 


Softcninj 


? (kpc) 




fv 


Name 




(M Q ) 


DM 


Gas 


DM 


Gas 






Al 


9171 x 2 


8 x 10 12 


7.85 x 10 8 


8.72 x 10 7 


5.39 


2.59 


0.1 





A2 


9171 x 2 


2 x 10 12 


1.96 x 10 8 


2.18 x 10 7 


3.40 


1.63 


0.1 





A3 


9171 x 2 


5 x 10 11 


4.91 x 10 7 


5.45 x 10 6 


2.14 


1.03 


0.1 





A4 


9171 x 2 


1 x 10 11 


9.81 x 10 6 


1.09 x 10 6 


1.25 


0.60 


0.1 





Bl 


9171 x 2 


8 x 10 12 


7.85 x 10 8 


8.72 x 10 7 


5.39 


2.59 


4 


0.9 


B2 


9171 x 2 


2 x 10 12 


1.96 x 10 8 


2.18 x 10 7 


3.40 


1.63 


4 


0.9 


B3 


9171 x 2 


5 x 10 11 


4.91 x 10 7 


5.45 x 10 6 


2.14 


1.03 


4 


0.9 


B4 


9171 x 2 


1 x 10 11 


9.81 x 10 6 


1.09 x 10 6 


1.25 


0.60 


4 


0.9 



Table 3 

Global photometric properties. 



Model 




V band 




K band 




A(B - R) 


Alog(Z/Z ) 


Name 


M B 


M v 


n 


r e 


M K 


n 


r e 


Alog(r) 


Alog(r) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Al 


-21.83 


-22.69 


4.18 


5.41 


-25.82 


3.96 


4.26 


-0.16 


-0.38 


A2 


-20.29 


-21.24 


3.17 


3.86 


-24.34 


3.06 


3.12 


-0.20 


-0.51 


A3 


-18.70 


-19.66 


3.41 


2.42 


-22.76 


3.39 


2.00 


-0.11 


-0.40 


A4 


-16.86 


-17.83 


3.51 


0.77 


-20.96 


3.42 


0.66 


-0.20 


-0.52 


Bl 


-21.63 


-22.58 


2.12 


5.45 


-25.82 


3.96 


4.26 


-0.33 


-0.63 


B2 


-20.19 


-21.14 


1.41 


4.32 


-24.24 


1.35 


2.86 


-0.17 


-0.26 


B3 


-18.66 


-19.59 


0.92 


3.63 


-22.64 


0.89 


3.45 


-0.10 


-0.15 


B4 


-16.72 


-17.63 


0.96 


1.45 


-20.56 


0.95 


1.40 


-0.08 


-0.11 



Note. — Col. (10): Luminosity weighted metallicity. 



D. Kawata 



Table 4 

Photometric properties and stellar populations within apetures. 



Model 




D <5 kpc 






D <99 kpc 






Name 


U - V 


V - K 


[Z/Hj 


Age 


U — V 


V - K 


[Z/Hj 


Age 


[Mg/Fe] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Al 


1.63 


3.31 


0.35 


11.3 


1.35 


3.10 


0.16 


11.6 


0.05 


A2 


1.58 


3.27 


0.33 


11.7 


1.35 


3.09 


0.14 


11.9 


-0.03 


A3 


1.51 


3.22 


0.27 


11.8 


1.34 


3.09 


0.13 


11.9 


-0.06 


A4 


1.48 


3.19 


0.20 


11.6 


1.34 


3.09 


0.15 


11.8 


-0.03 


Bl 


1.65 


3.32 


0.28 


11.9 


1.35 


3.10 


0.09 


12.0 


0.48 


B2 


1.54 


3.23 


0.15 


11.8 


1.32 


3.08 


0.03 


11.9 


0.48 


B3 


1.41 


3.11 


-0.04 


11.6 


1.29 


3.04 


-0.12 


11.7 


0.47 


B4 


1.26 


2.96 


-0.22 


11.4 


1.22 


3.01 


-0.25 


11.5 


0.47 



Note. — Col. (4) and (8): Luminosity weighted metallicity log(Z/Z Q ). Col. (5) and (9): 
Luminosity weighted age (Gyr). Col. (10): [Mg/Fe] in the r e ,s/2 aperture. 



Table 5 

SN rates. Units are SNu = SN(100yr)- 1 (10 10 L Q , s ). 



Model 


1.5 Gyr a 


2 Gyr b 


SNe la 


SNe II 


SNe la 


SNe II 


Al 


0.18 


0.47 


0.17 


0.74 


A2 


0.09 


0.38 


0.13 


0.37 


A3 


0.10 


0.43 


0.08 


0.47 


A4 


0.07 


0.23 


0.08 


0.31 


Bl 


0.00 


0.00 


0.00 


0.00 


B2 


0.00 


0.00 


0.00 


0.00 


B3 


0.00 


0.00 


0.00 


0.00 


B4 


0.00 


0.00 


0.00 


0.00 


E-S0 c 


0.18±0.06 


<0.02 







a W Avarage SNe rates for 1.5 (2) Gyr. 

c Observed SN rate in Caperrallo, Evans, & Tu- 
ratto (1999). 



